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Abstract 

s. 

, The Landau-Lifshitz fluctuating hydrodynamics is used to study the statis- 

Q , tical properties of the hnearized Kolmogorov flow. The relative simplicity of 

O ' 

this flow allows a detailed analysis of the fluctuation spectrum from near equi- 
^ \ librium regime up to the vicinity of the first convective instability threshold. 

It is shown that in the long time limit the flow behaves as an incompressible 
fluid, regardless of the value of the Reynolds number. This is not the case 
for the short time behavior where the incompressibility assumption leads in 
general to a wrong form of the static correlation functions, except near the 
instability threshold. The theoretical predictions are confirmed by numerical 
simulations of the full nonlinear fluctuating hydrodynamic equations. 
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I. INTRODUCTION 



A common theoretical approach for the study of fluctuations is the Landau-Lifshitz fluc- 
tuating hydrodynamics mainly because of its relative simplicity as compared to more 
fundamental approaches 0,^. Fluctuating hydrodynamics is a stochastic formulation of 
standard fluid mechanics. Spontaneous fluctuations of hydrodynamic variables are intro- 
duced into the transport equations by adding random components to the dissipative part 
of the pressure and heat fluxes. Since these fluxes are not conserved quantities, the cor- 
relations of the random terms are expected to be short-ranged and short-lived, so that on 
a hydrodynamic scale they are assumed to be Dirac-delta correlated. Their strengths are 
then chosen to yield the correct equilibrium thermodynamic fluctuations as derived from 
the Gibbs distribution. 

Fluctuating hydrodynamics has been used by various authors to study the statistical 
properties of simple fluids subjected to nonequilibrium constraints, such as temperature 
gradient or shear (for a review, see ref. 0). Recent light scattering results, 

obtained for systems under temperature gradient, have shown quantitative agreement with 
theoretical predictions @|. Quantitative agreements have also been demonstrated with re- 
sults based on particle simulations, both for systems under temperature gradient and 
shear |]TT|. A more important issue is obviously the role of fluctuations in the onset of hy- 
drodynamical instabilities, such as the convective instability arising in the Benard problem 

El. 

The macroscopic studies of sub-sonic hydrodynamical instabilities are usually based on 
the incompressibility assumption. However, as first pointed out by Zaitsev and Shliomis 
||13|| , this assumption is basically inconsistent with the very foundation of the fluctuating 
hydrodynamics formalism since it imposes fictitious correlations between the velocity com- 
ponents of the fluids. On the other hand, even in the absence of noise, the mathematical 
analysis of convective instabilities arising in compressible fluids proves to be quite involved 
|]1^,|T5|. One way to overcome this difficulty is to look for idealized models which, in spite 
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of their extreme simplicity, can nevertheless lead to hydrodynamical instabilities analog to 
those observed in real systems. Our main purpose in this article is to study the statistical 
properties of one such a model proposed some fifty years ago by Kolmogorov [0. As we 
will show, the periodic boundary conditions associated to this model allow detailed analysis 
of the fiuctuation spectrum from near equilibrium up to the vicinity of the first instability 
leading to convective rolls. 

The Kolmogorov fiow will be presented in the next section where some well known 
aspects of its macroscopic behavior are reviewed. The statistical properties of the model 
will be discussed in section III. We will show that the dynamic structure factor of the 
fiuid is practically not affected by nonequilibrium constraints. This is not the case for the 
velocity correlation functions which become long-ranged as soon as the system is driven out 
of equilibrium. Their amplitude is shown to diverges as one approaches to the convective 
instability threshold. Conclusions and perspectives will be presented in section IV. 



II. KOLMOGOROV FLOW 

Consider an isothermal flow in a rectangular box x Ly oriented along the main axes, 
that is{0<x<L^,0<?/<Lj/}. Periodic boundary conditions are assumed in both direc- 
tions and the flow is maintained through an external force fleld of the form 

Vext = Fq sai{2-nny/Ly)lx (1) 

where 1^ is the unit vector in the x direction. This model represents the so-called "Kol- 
mogorov flow" and it belongs to the wider class of two-dimensional negative effective eddy 



viscosity flows JT^. It is entirely characterized through the strength of the force fleld Fq, 
the parameter n, which controls the wave number of the forcing, and the aspect ratio ar, 
deflned as 

Or = Lx/Ly (2) 
In the following, we will mainly concentrate on the case n = 1. 



The hydro dynamic equations for this model read: 

If = -V . (pv) (3) 

^ -p(v- V)v - VP - V-o- + Fe,i (4) 

where p is the mass density, P the hydrostatic pressure and cr the two dimensional stress 
tensor: 



a. 



For simphcity, we shall assume that the shear and bulk viscosity coefficients, t] and (, are 
state independent, i.e. they are constant. It can then be easily checked that at the stationary 
state the pressure and the density are uniform in space (Pst = Pq, pst = Po) whereas the 
velocity profile is given by: 

= Mo sin (2 TT y/Ly) Ij; (6a) 

«o= ^ (6b) 

For small enough Fq, this stationary flow is stable. As we increase Fq, however, the flow may 
become unstable giving rise to rotating convective patterns. Other instabilities of increasing 
complexity may occur for larger values of Fq, culminating in a turbulent - like behavior 
||T8|^0[. In this paper we shall limit ourselves to the analysis of the system before its ffist 
instability. 

We still have to supply the momentum conservation equation, eq.(^, with an equation 
of state relating the pressure to the density (recall that the system is isothermal). In this 
section, we follow the traditional macroscopic analysis by assuming that the flow is incom- 
pressible, i.e. 

^ ^ _ du ^ dv _ ^ 
dx dy 

where u and v represent the x and y components of the velocity, respectively, i.e. v = 
ulx + vly. Relation (|^ implies a uniform density po throughout the system for all times, if 



initially so, as well as the existence of a scalar function ipi^j u)-, known as "stream function" , 
defined by the relations: 



u 



dip 
dy 



V = — 



dip 
dx 



(8) 



Scaling lengths by Ly, velocities by Uq and time by Ly/uo, the dimensionless equation for 
the stream function reads: 



dt dy dx dx dy 

where R represents the Reynolds number, 



Stt'^ R-^ cos{2 7iy) (9) 



R 



Uq Ly 



(10) 



(11) 



and V = 7]/ pis the kinematic viscosity. The stationary solution of (^) reads: 

ipst = - ^ cos (2 Try) 
Setting ip = ipst + ^V^, and linearizing (|^) around ipst, one gets 

= - si„(2.,)^<^ - 4.^ Bin (2.,)^ + R-^VHVH4.) (12) 
Owing to periodic boundary conditions, 6ip{x,y,t) can be expanded in Fourier series: 

oo 

5ijj{x,y,t)= exp{-2 piikyy) exp{-2nik^x/ar)5ijjk^^kyi't) , 

^i'k^,kAt) = dy exp {2 TT i ky y) — dx exp {2 7i i kj:x/ar) Sip{x, y, t) 
Jo ar Jo 

(13) 

Equation ([12D can then be transformed to 



dt 



-An'R-' {ki + ky')5ilJk^,ky 



+ TTk^ 



+ 27T^ ^"^^ 



^^j^hx 1 ky -\- 1 ^^^kx 5 ky 

5lPk,,ky + l + SlPk^^ky 



kl + ky^ 



(14) 
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where we have set 



In its general form, the analysis of this equation proves to be quite difficult On 
the other hand, if ipst is stable then, in the long time limit, the evolution of the system 
will be mainly governed by long wave length modes. Accordingly, we start our analysis by 
considering only the modes ky = , ±1, i.e. we assume that Sipk^,ky{t) ~ for \ky\ > 2 | |22[| . 
The equation (|Tip then reduces to a 3 x 3 matricial equation whose eigenvalues for kx ^ 
(the case k^ = is trivial) are found to be : 



Ai = -2rt^R'\l + 2kl) + T:pkl{l-kl)/{l + kl) + 47r2i?-2 



As = -2Ti^R-\l + 2kl) - Ti^2kl{l-kl)/{l + kl) + Ati^R-\ 

As = -At^^R'^{1 + kl). (16) 



It follows from ([iq ) that A2 and A3 are always negative, whereas there exists a critical 
value of the Reynolds number 

R,{k.x) = 2v^7r ] ^^L ; < < 1 (17) 
Vl - ~kl 

for which Ai becomes equal to zero, thus indicating the limit of the stability of the corre- 
sponding mode |2^. Clearly R^ is an increasing function of \kx\ so that the first modes to 
become unstable correspond to \kx\ = 1, provided the aspect ratio > 1. As — 1, 
Rc ^ 00, indicating that no instability can develop for perturbations of the same spatial 
periodicity as the applied force (see ref. [0]). In the following, we shall therefore concentrate 
on the case > 1. 

For ar = 2, relation (|1^ predicts a critical Reynolds number Rc = 57r-\/6/3 ~ 12.8255. 
Analytical calculations can still be handled when the modes ky = ±2, ±3 are taken into 
account as well, and lead to a critical Reynolds number Rc ~ 12.8736. The discrepancy thus 
remains smaller than 0.4%. Numerical evaluation of Rc performed with a total amount of 
103 modes shows no further discrepancy. We thus conclude that one can rely reasonably 



well on a "3-modes" approximation theory (that is SilJk^,ky(t) ~ for \ky\ > 2), as long as 
R is close to Re- We shall use this approximation in the next section to study the statistical 
properties of the system near its first instability. 



III. HYDRODYNAMIC FLUCTUATIONS 

To study the fluctuations spectrum, we first linearize the hydrodynamic equations ^ 
around the stationary state. Setting p = po + 6p, P = Pq + 6P and v = v^j + 5v, and 
following Landau and Lifshitz [jl|, the fluctuating hydrodynamic equations read: 
d6p (d5u d6v\ . ,^,d5p 

d 5v 

Po -j^ = -Po (vst • V) (5v - po (5v ■V)vst - V6P - V -6(7 (19) 

6(7 is the two dimensional fluctuating stress tensor: 

Sa,, = -v{j^+^-k,V. 5v) - C S,, V ■ 5v + S,, (20) 

where S is a random tensor whose elements {Sij} are Gaussian white noises with zero mean 
and covariances given by 

< ^.,,(r, t) SkAr', t') > = 2kBTo 5{t - t') S{r - r') [r/(<5g-<5f/ + ^f/^jj) + (C - vK'S^j] (21) 

We still have to specify the equation of state. Since the fluid is isothermal, we simply set 

6P = cl6p, (22) 

where Cg is the isothermal sound speed. Scaling lengths by Ly, time by Ly/cg, Sp by po and 
5v by the speed of sound, the dimensionless fluctuating equations in the Fourier space read: 

dSpk,,ky{t) 



dt 

dSuk^,ky{t) 



2TTi (k^dUk^^ky + ky^'^k^,ky) + R ^ P k^ , ky + 1 " ^Pk^,ky~l), (23) 

= ~TTeRi5Vk^,ky + l + SVk,,ky^l) + 7ieRh{SUk,,ky + l " 5Uk^,ky-l) 

- 4:iT'^e{kl + kl)5uk,,ky - 4:7r'^ a ekAkx5uk,,ky + ky5vk,,ky) 
+ 27iik^ 6pk^, ky + Fk^, ky {t) , (24) 
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- An"^ aeky{k^5uk,,ky + kySvk,,ky) + 27r i ky Spk,,ky + Gk,,kyit). (25) 
where a = (/rj and 

The functions Fk^^ky and Gk^,ky are Fourier components of the noise terms ; their covariances 
follow directly from eqs. (pO|, pT]): 



< Fk^,ky{t) Fk,,k'^{t') > =eA[{a + l)~kl + kl] 5^l^,^,5{t - t') , 

< Fk.., ky {t) Gk'^, k'y {t') > =eAahky 5^+^,0 - t') , 

< Gk^,kyit)Gk^^,k'y{t') > =eA[kl + {a + l)fcJfe,o'^(t - 0, (27) 



where k = {k^/ar, ky) and 

M = ttr po Ly'^ being the total mass of the system. 

The analysis of the above Langevin equations can be simplified somewhat by noticing 
that the quantity e must remain small if one wishes to remain within the limit of validity 
of the hydrodynamic regime Furthermore, as already mentioned in the introduction, 



in this article we limit ourselves to strictly sub-sonic flows, so that eR = Uq/cs << 1. We 
thus have at our disposal a natural small parameter which, however, have to be used with 
care since the solution of the Langevin equations - 125|) proves to be singular in the limit 
e — > 0. Moreover, it turns out that the behavior of the system is qualitatively independent 
of the value of the bulk viscosity coefficient so that, to avoid cumbersome notations, we 
simply set a = (recall that a = Q/rf). In any case, the calculations remain lengthy and 
tedious, so that here we concentrate mainly on the final results, giving only a brief sketch 
of the intermediate steps. 
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We pay particular attention to two quantities. First, the so-called scattering function, 
defined as the space-time Fourier transform of the density auto-correlation function: 

/+00 
dt exp (itot) < 5pkit) 5p_k(0) > (29) 
-oo 



< 5pk(t)5p_k(0) >= 1 1 drdr' exp{2mk- {r-r')} < 6p{r, t) 6p{r' , 0) > (30) 



where the integrals extend over the surface S" = x 1 of the system. Next, the space-time 
Fourier transform of the velocity auto-correlation function, defined in a similar fashion: 

/ + 00 
dt exp{iujt) < 5vk(t) ■ 5v_k(0) > (31) 
-oo 

as well as its static (equal time) counterpart: 

1 r+co 

< (5vk ■ 5v_k > = — / duj Wkiuj) (32) 

2 TT J-oo 

We found that, to dominant order in e, S]s_{u) is not affected by the nonequilibrium 
constraints, i.e. 

Sk{uj) = 2 1 + 0(e^R^) (33) 

where k"^ = {k^. + ky). We note that the scattering function exhibits only sound mode peaks 
(Brillouin lines). The absence of a purely dissipative mode around u ^ (Rayleigh line) is 
directly related to the fact that the Kolmogorov flow is "strictly" isothermal, i.e. there are 
no temperature (or entropy) fluctuations. On the other hand, the velocity auto-correlation 
function do exhibit a purely dissipative viscous regime around ~ 0, together with a sound 
regime located around u ^ ±27r|/c|. Here again we found that, to dominant order in e, the 
sound regime is not affected by the nonequilibrium constraints and behaves very much like 
the scattering function, eq. (p3|). One thus arrives at the conclusion that the nonequilibrium 
constraints affect mainly the behavior of the fluid near the origin ~ (the viscous regime). 

We first consider near equilibrium situations, limiting ourselves to relatively small val- 
ues of the Reynolds number R. In this case the Langevin eqs. (|2^ - |2^) can be solved 
perturbatively by expanding the hydrodynamic variables around the equilibrium. 



5p = 5peq + fiSpi + p'^ 5p2 + ■ . . (34a) 

(5v = 6\eq + /i 5vi + /i^ 5V2 + . . . (34b) 

where the subscript "eg" denotes equihbrium quantities and the parameter p is defined as 
p = R/2tt. After some tedious algebra, one gets for the static correlation function: 



< (5vk ■ 5v_k > — < ■ 5v_k > 



eq 



{R/2nf — ! , . l + 0(e^ (i?/2vr)2) (35) 

where < (5vk ■ 5v_k >eg = is the equilibrium contribution and, to simplify the presenta- 
tion, we have considered the case ky = 1. 

To check the validity of this result, we have solved numerically the Langevin eqs. (p3|- ^5|). 
The traditional procedure consists of simulating the corresponding stochastic processes and 
using the hydrodynamical sample paths (time series) so obtained to construct the various 
correlation functions. Unavoidable for nonlinear problems, this procedure is quite simple to 
set up but requires very long runs in order to get reliable statistics. Alternatively, one can 
solve directly the equations governing the evolution of the correlation functions which can 



be obtained easily from the underlying Langevin equations |26|. This latter procedure is 
accurate (no need for statistics) and quite fast but it is, of course, limited to linear problems. 
We have used both techniques, the former to simulate the full nonlinear hydrodynamic 
equations (^ ^) with noisy stress tensor, and the latter to study the statistical properties of 
the (linear) Langevin equations (P3| - limiting ourselves to the first 41 ky modes (that 
is Spk.,ky{t) = S^k^,ky{t) ^ for I A;^^ I > 21). 

In figure 1 we have presented the static velocity auto-correlation function, as given by 
the relation (|35|) , together with the corresponding numerical solution. As it can be seen. 



quantitative agreement is demonstrated for R < 4 but discrepancies gradually appear as we 
consider larger values of the Reynolds number. This is to be expected since the validity of 
the relation (^) can only be guaranteed for "small" values of the Reynolds number. 
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Before discussing the behavior of the system for larger values of R {< Rc) it is instructive 
to study the properties of the static correlation function in real space, < 6\{r) ■ (5v(0) >. 
This can be obtained by summing the product < (5vk ■ 5v_k > x exp[27rz(x + yky)] over 
{kx, ky). Analytical calculations, however, prove to be extremely difficult to handle for the 
general case. We therefore limit ourselves to a special case where only one of the wavenumber 
is summed over, the other being held fixed. Specifically, we simply set A;^ = to obtain: 

< 5v,, ■ 5v_,, > -2A = {R/2nf ^''^ [l + 0(e^ {R/27,f)\ (36) 

Note that setting ky = Q is equivalent to taking the spatial average over the y direction (cfr. 
eq. (|30|)), so that relation (^) holds only for kx 7^ 0. In fact, 5vo, o(i) = since the hnear 
momentum of the center of mass is a conserved quantity. With this restriction in mind, the 
summation can be performed in a straightforward manner to give 



AR^ ttr 

Stt 



< 5v(x) • 5v(0) > -2A[6{x) - ar 
cosh[v^ 7r (x — ar/2)] cosh[2 n (x — ar/2)] 



(37) 



sinh(7r 0,71/2) sinh(7rar) Trar 

where the second term on the left hand side is the equilibrium contribution |^. Note the 
presence of a constant term in both equilibrium and nonequilibrium (right hand side) part 
which ensures the conservation of the linear momentum. 

The nonequilibrium contribution to < 6v{x) ■ 5v(0) > exhibits long-range correlations 
since the correlation length is clearly of the order of system's size. This is shown in figure 
(2) for i? = 3 where quantitative agreement with numerical results is observed. The exis- 
tence of long-range correlations is generic for fiuids under shear constraints and have been 



predicted by several authors [0,0, and confirmed by both microscopic |jTT| and lattice-gas 
automata simulations On the other hand, experimental evidence has been so far re- 
ported only for fiuids under a temperature gradient, where quantitative agreement with 
fiuctuating hydrodynamics was demonstrated p]. 

Let us now consider far from equilibrium situations. As pointed out in the last section, 
for R close to Rc one can reasonably well rely on the " 3- modes" approximation theory (that 
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is Spk^^ky{t) = Svk^^ky{t) ~ for \ky\ > 2). As a consequence, equations - 0) reduce to 
a system of nine coupled linear Langevin equations for each fixed which, for consistency, 
must be limited to \kx\ < a^. The calculations can nevertheless be done, leading to the 
following expression for the static velocity auto-correlation function: 



< 5vw ■ (5v_w > -2A 



1 + Oie'W) (38) 

2{Rl-R^)[l + 2kl] 

where the second term on the left hand side is the equilibrium contribution, Rc{kx) is given 
by eq. (\[7\} and ky = 1. The nonequilibrium part diverges as i? — ^ Rdk^ = cir), but then, of 
course, the linearized Langevin equations ceases to be valid. In figure 3 we have represented 
the result (|38D for increasing values of R, together with the numerical solution of the linear 
Langevin equations (|2^- ^) as well as the results obtained by simulation of the full nonlinear 
hydrodynamic equations (^, |^) with noisy stress tensor. Quantitative agreement is observed 
for values of R up to 12, but significant discrepancies start to show up as i? ^ -Rc(~ 12.87) 



where the linearized theory leads to diverging correlation function (cfr. eq. (|38|)). This is 
not the case for the correlation function based on the full nonlinear equations which seems 
to exhibit a maximum around Rc. It should however be noted that, due to slowing down 
of the relaxation of the " critical" Fourier modes, statistical errors are quite important for R 
close to Rc (about 15% for the last four data), so that no definitive conclusion can be drawn 
at this stage. In any case, the analysis of the statistical properties of the nonlinear regime 
is beyond the scope of the present work. 

IV. VALIDITY OF THE INCOMPRESSIBILITY ASSUMPTION 

The macroscopic studies of sub-sonic hydrodynamical instabilities are based on the in- 
compressibility assumption which is fundamentally inconsistent with the very foundation 
of the fluctuating hydrodynamics formalism [|1^]. For instance, it is easy to show that at 
equilibrium {R = 0), one has: 

UUu^U _ {^'-^-'kl)\ 16.^ klkl 



^k(w)eg (lu^ -A + 16 Tc^klk? 



1 + 0{e') (39) 
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where Uk^u) and Vk(co') are the space-time Fourier transforms of < Su{T,t)Su{r',0) > and 
<6v{r,t)6v{r',0) > respectively. On the other hand, the incompressibihty assumption, eq. 
(0), imphes 

Ui,ico)rnc/V^iu)inc = kl/kl (40) 

where the subscript " inc" refers to incompressible fluids. Except near the origin uj ^ (long 
time limit), this result is clearly in contradiction with the correct equilibrium form, eq. (p9|). 
In particular, the equilibrium static auto-correlations are independent of the wave vector, 

< 5nk 5M-k >eq = < SVk SV-y, >eq = A (41) 

whereas relation ( ^OD leads to < ^Wk^w^k >mc / < '^'^k'^'W-k >inc= ky / k^. 

The situation is somewhat different for the nonequilibrium case. As we have shown in 
the previous section, to dominant order in e both the scattering function and the sound 
regime of the velocity auto-correlation function assume their equilibrium form, regardless 
of the value of the Reynolds number. The nonequilibrium constraints thus affect mainly 
the behavior of the fluid near the origin u ^ (the viscous regime). This result has an 
interesting consequence. It suggests that, as far as the nonequilibrium properties of the 
fluid are concerned, one may rely on the incompressibihty assumption, eq. (^, since the 
compressibility of the fluid affects mainly the sound modes which are well separated from 
the purely viscous modes (recall that e is small). Analytical calculations confirm the above 
arguments and lead to the following relation: 

^k(^)-j7k(^)e. ^gr .^.1 (42) 

Vic(^) - V^k(^)eg fc2 L ^ ^ ^ 

where both the numerator and the denominator in the left hand side prove to assume a 
Lorentzian shape, sharply peaked around the origin (the width is of the order of e), since 
the sound regime cancels out. Nevertheless, because of the presence of the equilibrium 
contributions, the relation (03) is still in contradiction with the incompressibihty condition. 



eq. (HO). There exist, however, two different situations where this objection can be ruled 
out. 
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First, near the origin {u ~ 0) where the fluid satisfies the incompressibihty condition 
already at equilibrium, i.e. Uk{Lj)eqkl ~ V]^{uj)eqky. Obviously, this situation concerns only 
the long time behavior of the fiuid. For instance, the static correlation functions obey 



y 



l + 0{e'R') (43) 



< Svk6v-k > - A kl 
which contradicts the incompressibihty condition, eq. (0). 

A more interesting situation concerns the behavior of the fiuid near the instability where 
it can be shown that in the limit R ^ Rc both the static and dynamic velocity correlation 
functions behave as 0{{Rl — R'^)~^). In other words, for R "close" enough to Rc, equilibrium 
contributions become negligible so that the fiuid behaves basically in an incompressible way. 

It should however be realized that this appealing conclusion is based on the linearized 



Langevin equations (|2^ - ^) which are not valid near the convective instability. The study 
of the statistical properties of the system in the critical regime requires a nonlinear analysis 
of the fiuctuating equations which is beyond the scope of the present work and will be 
reported elsewhere. Instead, we resort here to numerical analysis only. More precisely, we 
have simulated the full nonlinear fiuctuating hydrodynamic equations to obtain the ratio of 
the X and y components of the static velocity auto-correlation function for several values of 
the Reynolds number. The results are depicted in figure 4 for kx = ky = 1, ar = 2, so that 
the expected value of this ratio for an incompressible fiuid is 4. This is precisely what we 
observe, but only for values of -R > 12.8 (recall that Rc ~ 12.87), a domain which is beyond 
the validity of linearized hydrodynamic equations (cfr. previous section). 



V. CONCLUDING REMARKS 

Our main purpose in this article was the study of the statistical properties of the lin- 
earized Kolmogorov fiow. The simplicity of this model allows a detailed analysis of the 
fiuctuations spectrum from near equilibrium up to the vicinity of the first instability lead- 
ing to convective rolls. For this latter case, the analytical calculations were based on a 
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"3-modes" approximation theory, which consists in retaining only the Fourier modes with 
wavenumber \ky\ < 1, while for the former case we have set up a perturbation scheme 
around the equilibrium. Extensive numerical calculations allow us to emphasize clearly the 
limit of validity of both regimes. In particular, we have shown that the " 3-modes" approx- 
imation theory holds already for i? > 8 and leads to a diverging velocity auto-correlation 
function as R ^ Re- On the other hand, the simulation of the full nonlinear fluctuating 
hydrodynamic equations indicates that the validity of linearized hydrodynamic equations 
can be guaranteed for Reynolds number as high as 12 {Rc ~ 12.87). 

Another interesting result concerns the validity of the incompressibility assumption which 
greatly simplifies the mathematical analysis of the problem. The compressibility of a fluid 
affects mostly fast sound modes whereas the dynamic of the system near an instability is 
expected to be governed mainly by dissipative slow modes. This intuitive argument has 
been used by many authors who have considered fluctuating incompressible hydrodynamic 
equations, or even directly the corresponding normal form amplitude equations to which 
they added random noise terms In these approaches, the characteristics of the noise 



terms cannot be related to equilibrium statistical properties of the fluid and thus remain 
arbitrary. A more satisfactory approach would be to start with the compressible fluctuating 
hydrodynamic equations. Such a procedure, however, proves to be extremely difficult mainly 
because of the boundary value problem. To our knowledge, the only attempt in this direction 



has been made by Schmitz and Cohen for the case of Benard instability [Tj| . Concentrating 
on the behavior of a small layer in the bulk, these authors have succeeded in deriving the 
linearized fluctuating equations close to the convective instability. 

Here again, the relative simplicity of the Kolmogorov flow allows some further progress 
in this important issue. In this respect, we have shown that in the long time limit the 
flow behaves as an incompressible fluid, regardless of the value of the Reynolds number. 
This, however, proves to be not the case for the short time behavior. In particular, the 
incompressibility assumption leads in general to a wrong form of the static correlation func- 
tions. The only exception is near the convective instability, where we have shown that the 

15 



incompressibility assumption remains valid. 

The problem with this conclusion is that the linearized fluctuating hydrodynamic equa- 
tions, on which this whole paper is based, are no more valid close to the instability threshold. 
Although extensive numerical simulations have basically conflrmed our predictions, a full an- 
swer to this important problem requires nevertheless a nonlinear analysis of the fluctuating 
Kolmogorov flow. Work on this direction is in progress. 
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FIGURES 

FIG. 1. Fourier transform of the nonequilibrium part of the static velocity auto-correlation 

function, normalized by tlic corresponding equilibrium part, as a function of the Reynolds number. 
The solid curve represents the theoretical prediction, as given by eq. (35), whereas the diamonds 
correspond to results obtained by numerical simulation of the linear Langevin equations (23 - 25). 
The parameters are : e = 10"^, A = 10^^ (defined by eq. (28)), ky = l,kx = 1 (i.e. kx = 1/2). 
The estimated statistical errors are less than 4%. 

FIG. 2. Nonequilibrium part of the static velocity auto-correlation function, normalized by 
the parameter A (defined by eq. (28)), as a function of the spatial coordinate x. The solid curve 
represents the theoretical prediction, as given by eq. (37), whereas the diamonds correspond to 
results obtained by numerical simulation of the linear Langevin equations (23 - 25) . The Reynolds 
number is set to i? = 3 and the other parameters are as given in the caption of figure (1). 

FIG. 3. Fourier transform of the nonequilibrium part of the static velocity auto-correlation 
function, normalized by the corresponding equilibrium part, as a function of the Reynolds number. 
The solid curve represents the prediction based on the " 3- modes" approximation theory, eq. (38) , 
whereas crosses and diamonds correspond to numerical results obtained respectively by the simu- 
lation of linear and nonlinear Langevin equations, eqs. (3, 4). The parameters are as given in the 
caption of figure (1). The estimated statistical errors are about 15% for the last four data points. 

FIG. 4. Ratio of x and y component of the velocity auto-correlation function in Fourier space, 
as a function of the Reynolds number. The parameters are : e = 10~^,A = 10~^ (defined by 
eq. (28)), ky = l,kx = I (i.e. kx = 1/2). The dashed line represents the expected ratio for an 
incompressible fluid, whereas the diamonds correspond to results obtained by the simulation of the 
full nonlinear Langevin equations, eq. (3, 4). The estimated statistical errors do not exceed 8%. 
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